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Abstract 



In this letter we present a computation of the phase induced by test masses 
of different geometry, in the framework of non-newtonian gravitation, on an 
ideal separated arms atom interferometer. We deduce the related limits on 
the non-newtonian gravitational strength in the sub-millimeter region for the 
potential range. These limits would be comparable with the best existing 
experimental limits but with the advantage of using a microscopic probe. 



In this first section we will recall briefly the main lines of the theoretical motivations of 
this kind of calculations and their interest and the present state-of-the-art of experiments in 
this field will be summarized. For a more detailed description of both theoretical and exper- 
imental point of view the reader is referred to [|IJ and 0. The basics of atom interferometry 
will be then presented and, after some geometrical considerations, the signal for different 
kinds of potential will be obtained. Finally, we will derive the limits on the non-newtonian 
gravitational strength that could be obtained by an ideal experiment and we will compare 
the results to present limits. 

In the attempt to unify the description of all known forces, two types of extensions of 
newtonian gravity are usually made [|]]. First, one can postulate that a new force mediated 
by massive scalar bosons exists and that it gives rise to a Yukawa potential so that the 
total gravitational potential between point-like particles of mass m\ and rri2 separated by a 
distance r reads: 



In the preceding expression G denotes the gravitational constant and the first term is the 
usual newtonian potential. The correction term is thus characterized by its dimensionless 
strength ay and its range Ay which is essentially inversely proportional to the boson mass. 
A peculiar feature of such Yukawa potential is that the exponential factor saturates when 
the mutual distance r tends to zero: if its strength ay is small compared to 1 then the 
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Yukawa correction remains always a small correction to newtonian gravity which is itself a 
very weak interaction compared to the other known forces. 

On the other hand, some recent theoretical developments account for such weakness 
supposing that there exists n extra dimensions in which only gravity propagates so that 
the main part of it is lost for the usual 3D-space [[| . The extra dimensions are compact so 
that their existence can be felt only to mutual distances smaller to a compactification radius 
A n : for distances much greater than the compactification radius then the reminder of the 
gravitation interaction in the extra dimensions manifests itself as a Yukawa-like potential 
with a range of the order of the compactification radius A n and a strength a n of a few 
units. It thus quickly decreases leaving just the weak long range usual inverse square law 
interaction. On the contrary, at short mutual distances, it follows from Gauss theorem in 
3 + n spatial dimensions that the gravitational interaction V n (r) behaves as r~( 1+n ) so that 
it conveniently reads: 

V n (r) = -^pa n ^J forr«X n (2) 

It should be noticed that this potential does not saturates as the Yukawa potential at 
zero distances and "produces dramatic deviation from newtonian gravity" [l]]. The model 
predict a compactification radius A2 ~ 10~ 4 — 10~ 3 m which make it a matter of an 
experimental test. 

Up to now, there is experimental evidence that the inverse square law for gravitational 
interactions is well tested for separations r between 10 _2 m and 10 15 m where a is known 
to be smaller than 10~ 4 . A detailed description of the experimental works can be found in 
reference . The best limits on a in the region around A ~ 10~ 4 m have been very recently 
provided by the experiment described in reference ||. This experiment relies on a torque 
pendulum, the motion of which would be perturbed in a specific way by non-newtonian 
potentials from a test mass. Separation between 200/im and 10mm have been achieved with 
a null result. The present day constraint, therefore, for A ~ 10 m is a < 100 0]. 

Up to now all the performed experiments share common features from technical con- 
straints and scaling considerations. In particular, they are macroscopic in the sense that 
both the probe and test masses dimensions are in the centimeter range. 

In the following, we will show that atom interferometry, that uses a microscopic probe, 
could explore different geometries and distances ranges. So let us now turn to the basics of 
atom interferometry. 



II. ATOM INTERFEROMETRY 

Atom interferometry aims at making interferences with de Broglie waves associated to 
the external motion of massive particles. It was born in the late 80's when mechanical 
effects of light were extensively studied which lead the principal contributor, W. D. Philips, 
S. Chu and C. Cohen Tannoudji to Nobel price in 1997. The main difference with optical 
interferometry lies in the fact that the particle used (atoms or molecules) possess internal 
degree of freedom and mass. That make the interferometer sensitive to external fields and 
inertial effects we are concerned with. For an introductory review the reader is referred 
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to IJ and for more details to 0. As we want to probe distance dependent potential we 
will concentrate on an idealized separated arms interferometer as outlined in Fig. [l|. An 
atom beam is emitted from a source S and passes through a three identical beam splitters 
G\, 2, 3 separated by a distance L. At each beam splitter, each incoming beam is coherently 
divided in two parts separated by an angle 9 (labels A, O, O' , B on Fig. [I]). Only one of 
the closed paths is shown on Fig. p] at the end of which the interference pattern is recorded 
on a detector D. The detection scheme obviously depends on the particles used such as 
alkaline (Li, Na, Kb, Cs) atoms or molecules, earth-alkaline metastable atoms (Mg*, Ca*), 
metastable rare gas atoms (He*, Ar*, Ne*), molecules (I 2 , Na 2 , C 60 ). 

In the same way, the actual shape of the interferometer depends on the very technique 
used to construct the beam splitters. The example chosen above corresponds to grating 
interferometers: atoms, whose de Broglie wavelength is Xjb, interacting with a modulated 
perturbation of period A, are diffracted at a typical diffraction angle 9 ~ ^4 a . Such a 
perturbation can be for example a material grating or a laser standing wave. 

In ideal conditions, the signal recorded by the detector is then I(5<f)) = Iq cos 5<p where 
Iq equals to the atom flux times the interaction time and 5(f) is the phase difference between 
the two arms. 

In the following, to evaluate the phase difference induced by gravitational interactions, 
we will consider a semi-classical situation in which the atom velocity v and mass M are high 
enough so that its de Broglie wavelength X^b = ^ is small compared to any relevant length 
characterizing the potential V(t, t) existing in between the beam splitters. This potential 
will also be considered as a perturbation of the free motion of the particles which thus 
fly in straight line between the source, beam splitters and detector. Then, the phase shift 
accumulated by the particle with respect to the unperturbed case can be written as [§], 0: 

where r c / is the classical unperturbed path of the particle. The phase difference simply 
reads: 

50=1 j ^V(r,t(r))-i J ±V(T,t(r))= f ^V(v,t(v)) (4) 

ao'b aob ao'boa 

To effectively compute this phase difference and thus the expected signal one has specify 
the geometry of the test mass to evaluate the generated gravitation potential. 



III. GEOMETRY CONSIDERATIONS AND PHASE CALCULATIONS 

As explained in reference [0, one should use small size test mass to get as near as possible 
and thus increase the relative sensitivity to non-newtonian potentials. Nevertheless, as we 
investigate a totally different experimental technique, new possibilities are opened. Obvi- 
ously, a linear (ID) test mass parallel to the atomic trajectory is preferable to a point-like 
(OD) one as it increases the interacting time with the probe. One the other hand the use 
of a plane (2D) test mass is unfavourable for two reasons pi. First, as the potentials of 



interest are of short range, only a stripe of width ~ A parallel to the atomic trajectory will 
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contribute significantly to the interactions so that the ID case give a reasonable approxi- 
mation. Second, an atom near a surface is subjected to the Van der Waals potential that 
will screen the other ones. On the contrary, one might expect that the electrostatic image 
of the atom by a ID-wire is weaker than by a 2D-plane leading to a smaller Van der Waals 
interaction that thus will be neglected in the following. To be exhaustive, one should go 
into the higher dimensions cases. The 3D case of typical size R is also unfavourable because 
only a volume RX 2 has to be taken into account. However, as suggested in Eq. [| the use 
of time-dependent potentials (e.g. an oscillating test mass) open the way to fruitful phase 
sensitive detection techniques. 

We thus chose for the calculations an idealized wire test mass of linear density yu, of 
length L, located between the two first beam splitters at a distance d from the lower partial 
beam and parallel to it (see Fig. |1|). Let us take G2 as the x-axis of a reference frame and 
its intersection with the z-axis as the origin. Assuming a Yukawa extra interaction, the 
potential at the position (x, z) is: 

a Y GM , . 1 



V Y D (x,z) = . 1 fi du exp(- — J(d + x) 2 + u 2 ) (5) 

J-L-z ^(d + x) 2 + U 2 Ay V 

In the preceding expression, |<Z+x| is of the order of Ay, thus in the sub-millimeter range, 
whereas L is the length of the test mass which is related to the length of the apparatus 
which is commonly in the meter range. So, except for negligible fringe effects, the limits 
of the integral may be rejected to infinity as the exponential term to be integrated quickly 
decreases. Then, with the variable changes v = and v = sinhw: 

/+00 \jj _|_ nA \fj _|_ T | 

dw exp(-^^coshw) = —ay GMfi 2i^o(-^— r- — -) (6) 
-00 Ay Ay 

where Kq is the modified Bessel function of first kind . 

Between the two first gratings (z < 0), the two partial beams can be parameterized as 
x 2 (z) = and X\(z) = e(l + r)- So the phase delay accumulated along the lower path 
simply reads fa = G ^ v m 2cxyKo(-^), where m = Lfi is simply the test mass. The phase 
delay accumulated on the upper path (pi contains the mean value of the K function which 
has a rather complicated analytical expression involving hypergeometric and other Bessel 
functions. For sake of simplicity, we will consider only the limiting case where the beam 
separation e is much greater than the interaction range A. Then the potential is practically 
negligible everywhere along the upper path and so for <fii. The phase difference is then 
simply 50 = 2 — 0: 



GMm „ , d 



2a Y K (—) if e > Ay (7) 



hv A 



Y 



^From Ref. |TT| we get the following expansions: <50^? ~ yj je a as ^ <C 1 and S^if ~ In ^ 



as 4 > 1 



The same kind of calculations can be done for the case of extra dimensions. Some 
difficulties arise because the law we have (Eq. |^) is valid only for distances shorter than the 
compactification radius A n . For longer distances, the potential decreases exponentially fast 
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in a Yukawa fashion j|] and should give a negligible contribution. Let's make the assumption 
that d <C A n . As done before we will suppose that e ^> A to neglect the contribution of the 
upper path. We will get an estimate of the effect on the lower path restricting the integration 
range to a distance A n around the position of the particle. Under these assumptions, and 
neglecting fringe effects, the potential arising from n > extra dimensions reads: 

rriD, \ + f" a n GM X n 
v n 0, z ) = / r. V du 



^/(d + x) 2 + u 2 .^(d + xY + u 2 



a n GMfi(—^-) n [ dv{l + v 2 ) 
\d + x\ J 



l + n 
2 



\d+x\ 



where 2 -Pi denotes the hypergeometric function |13| . 
The phase difference then reads: 

GMm .Ans.ii „,1 1 + n 3 ,A n , 9 , „ 
™ = ImT n ( T } 2 l( 2 ' ~2~' 2 ' } ?/ 6 >> n andn>0 ( 9 ) 

and, in particular, 



One can show that, assuming that m is a constant, in the region d < A were the model 
holds that ctyjf ~ (^)" when d < A. 

Nevertheless, the newtonian potential n = has infinite range and thus must be treated 
separately. As the distance between the particle and the wire is much smaller than its 
length, one can take the infinite length of the wire limit [nj]. So, from Gauss theorem, the 
gravitational field is g(r) = —2 Gfi ^ = —V (2 Gfi In ^ so that the potential is null on the 
lower path. The potential energy is then: 

V^ wton {x) = 2GM^{\ + ^) (11) 

The newtonian phase difference is then: 



GMm „ , ,„ e 

'newton = ~ ^ — v- ■ rf 



1D 21n(1 + r) (12) 



With all these quantitative expressions for the different phases we can now turn to 
numerical estimates for a reference atom interferometer. 
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IV. NUMERICAL RESULTS 



In all cases, the phase difference is G ^ I v m times a numerical factor which can be inferred 
from dimensional considerations. The numerical factor depends on the particular geometry 
used and the presumed force law. It should be noticed that this numerical factor involves 
dimensionless ratios as a, 4 that can be great. Anyway, let us first give a numerical estimate 
of the pre-factor, taking a slowed cesium beam v = 10 m.s _1 ,M = 2 10~ 25 kg, and a wire 
lm x lOOyum x 100/im of density equal to 20 (gold) so that m = 2 10~ 4 kg. The pre-factor 
then amounts to some 10~ 6 . The detection limit 5<f) m i n = 10~ 3 rad is commonly accepted 
for atom interferometry (see different contributions in [7j]). It first confirms the known fact 
that the metrology on the regular newtonian gravity (i.e. a measurement of G) is practically 
impossible as even with d = 10//m, ln(l + §) will amount to some units so ^newton (Eq-|T2"D 
will be to small to be detected. 

One can anyway evaluate the upper bounds on the strength of the non-newtonian inter- 
action in the case of a hypothetical null experiment based on our calculation. The results are 
shown in Fig. |2] together with limits taken from Ref . |4[] . They were obtained assuming that 
the actual wire diameter equals the distance from the wire to the atoms d so that the ID 
model is at the limit of acceptability. The upper part corresponds to d = 100/im so that, for 
A < d, the wire is beyond the compactification radius for the models with extra dimensions. 
The curves have then been extrapolated the curve given by the Yukawa potential, d = 10/im 
is shown in the lower part for which such a problem do not arise since we have chosen the 
same parameters range as in Ref. 0j for sake of comparison. In the same way, the beam 
separation e ~ 10 _3 m with the chosen parameter set. The limits are thus overestimated in 
the somehow low interest region A > 10~ 3 m. 

While a real experiment based on this computation looks very difficult for several reason, 
the figure clearly shows that there is an interest in further studies of the methodology we 
described. For short ranges, A < 10 _4 m, the method appear particularly efficient even if one 
uses a less sophisticated atomic source such as a supersonic beam for which v ~ 10 3 m.s -1 . 
In the end, the main free parameter in the simulations is d and has to be chosen with 
respect to the potential model of interest and its suspected range. When d decreases the 
different potentials increase but the test mass scales typically as d 2 and vice versa. For 
the Yukawa case, an optimum is numerically, and not surprisingly, found for d ~ A for the 
typical parameters ranges we are concerned with. On the contrary for n extra dimensions 
the overall phase scales as d 2 (^j ~ d 2 ~ n . The case n = 2 thus exhibits a weak dependency 
on the choice of the wire diameter whereas one should use as small as possible a wire if 
n > 3. On the contrary, if the n — 1 case is under investigation, one should use the largest 
acceptable wire that is d ~ A. These behaviors can be checked on the figure where the limits 
for n — 1, 2, 3 get worse, remain roughly unchanged or improve from d = 100/xm (upper 
graph) to d — 10 fj/m (lower graph). 

V. CONCLUSIONS 

The results presented here concerning non-newtonian gravitation are based on a com- 
pletely different method than the existing ones that use macroscopic probes. Our computa- 
tion show that atom interferometry could provide limits on the strength a of several extra 
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potentials comparable with the existing ones in a reasonable integration time, especially if 
they are of short range A < 10~ 4 m. More, it can be easily extended to test a composition 
dependent interaction ("fifth-force") using, for example, as test masses two wires of different 
materials or two isotopes such as 6 Li and 7 Li for the probe atom. Using polarized beams 
spin dependent potentials could be also taken into account. 

A detailed study of an actual experimental set-up is far beyond the scope of such a letter; 
a rapid realization of such an experiment looks anyway difficult because of some experimental 
critical points. In particular the short range interaction implies high collimated beams which 
results in low count rates and moreover to increase sensitivity challenging beam slowing 
techniques are necessary. 

We are anyway confident that, being atom optics a fast developing field, in the near 
future most of the experimental difficult points will be clarified and atom interferometry 
will contribute to the tests of non-newtonian gravitation. 
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FIGURES 



FIG. 1. General scheme of a separated arms interferometer of overall length 2 x L and 
maximum separation e. S and D represent respectively the source and the detector. G\, G2 and 
G3 are the beam splitters that divide and recombine the atomic beam. The wire, of diameter d 
is set parallel and at the distance d of the first part of the lower arm. As an illustrative example 
pp. 1-83, a three gratings Mach-Zehnder interferometer of D. Pritchard's group, MIT, uses 
nanofabricated gratings of period A = 200nm separated by L = 0.6m. The source is a rare gaz 
seeded supersonic beam of sodium. Their de Broglie wavelength A^_b = ttj ~ 16pm. It corresponds 
to a diffraction angle 9 = A^s/A about 80/irad. The beam spacing e ~ 50fj,m is big enough so 
that a septum can be inserted between the two arms. The signal is then recorded on an hot wire 
detector. 

FIG. 2. Comparison of results presented in Ref. [|j (heavy lines) and detection limits of an 
atom interferometer assuming a Yukawa potential (thin line) or n extra dimensions n = 1 (dash), 
n = 2 (dot) or n = 3 (dot dash). The vertical axis, generically labeled \a\, represents either ay 
either 01,2,3 depending on the chosen scenario. Results in Q assume a Yukawa extra potential. The 
upper and lower parts correspond respectively to different wire diameters and beam/wire distances 
d = 100/xm and d = 10/um. In the former case, the results for extra dimensions models have been 
extrapolated to the Yukawa model by a dotted line in the region d > A. 
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